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We evaluate the two point functions of the stress energy from the largest string simulations 
carried out so far. The two point functions are used to calculate the cosmic microwave background 
(CMB) and cold dark matter (CDM) power spectra from local cosmic string models for structure 
formation. We find that our spectra differ significantly from those previously calculated for both 
global and local defects. We find a higher Doppler peak at I = 400 — 600 and a less severe 
bias problem than for global defects. Spectra were obtained for a variety of network energy-decay 
mechanisms. 
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1. INTRODUCTION 

The next few years will see a great increase in the accuracy of the data mapping 
the CMB temperature and large scale structure (LSS) due to satellite, balloon and 
ground based experiments. Similarly accurate theoretical predictions from inflation 
and defect models will be required for direct comparison with the new data. While 
calculations of structure formation from inflationary models are relatively straight- 
forward, defects, being highly non- linear objects and active sources of perturbations, 
pose a more challenging problem. 

To this end we report on what we believe to be the most complete numerical 
treatment of local cosmic strings to date. The method we use is based on that 
described in |l]] where it was shown how to use the two point functions (known 
as the unequal time correlators, or UETCs) measured form defect simulations to 
calculate the CMB and LSS power spectra. The method was applied to global 
defects and the results were discouraging for defect models. 

In this work we aimed to verify these results in local theories where the calcula- 
tions are complicated somewhat by considerations on energy conservation and have 
only been permitted following recent advances in the computer technology at our 
disposal. The method we used is considerably different from that used in previous 
local string calculations based on simulations and we believe that the UETC 
method is justified by causality and scaling arguments in dramatically extending 
the dynamical range of defect simulations. We also find that some of the our UETCs 
are approximated quite well by those based on an analytical model Q but others 
differ quite significantly and in j{| we |>how how this model can be elaborated to 
explain the features observed in the simulations. 



2. THE CORRELATORS 

The perturbations arising form scalar, vector and tensor (SVT) modes are com- 
pletely decoupled and evolve separately. To calculate the CMB and LSS power 
spectra efficiently, it is therefore convenient to decompose the source of the pertur- 
bations into irreducible components under rotations. 

Let M „(x) be the defect stress-energy tensor. We may Fourier analyze it 

e„„(x) - J d 3 fcO^(k)e lk x (1) 

and decompose its Fourier components as: 

Ooo = P d (2) 
OJ = ikv d + uf (3) 

&ij = y d k 3 + (kk - ^Si^j n s + 

i^kuJ + kjUY) + n£ (4) 

with Uujf = 0, kmY = 0, kmf 3 = 0, and Hp = 0. The variables {p d , v d lP d , n s } 
are the scalars, the vectors, and IT^ the tensors. 

The unequal time correlators of these components will completely specify the 
evolution of the network's stress-energy. All correlators between modes at (k, rj) 
and (k',77') will be proportional to S(k. — k') due to translational invariance. We 
shall drop this factor in all formulae. The correlators can also be functions of k 
alone, due to isotropy. Since conjugation corresponds to k — > — k isotropy implies 
that the correlators must be real. Because of incoherence the correlators will be 
generic functions of r/ and 77'. 

On top of this, the form of the S+V+T decomposition fixes further the form of 
the correlators. One can always write down the most general form of a correlator, 
and then contract the result with ki or 8{j, wherever appropriate, to obtain further 
conditions. Proceeding in this way we can show that cross correlators involving 
components of different type (S, V, or T) must be zero. Furthermore one has that 
the only non-zero correlators are the 10 scalar correlator functions 
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and the single tensor correlator function / . 

In general these functions are functions of (/«, 77,77'), and this is indeed the case 
during the matter radiation transition. However well into the matter and radiation 
epochs there is scaling, and these functions may be written as: 

f XY (k, V ,r 1 ') = FXY ^l X ' ) (7) 

v w 

where XY represents any pair of superscripts considered above, and x = kr] and 
x' = krf . The purpose of this work is the measurement of the 14 functions of two 
variables F XY (x, x'). 



3. NUMERICAL DETERMINATION OF THE CORRELATORS 

Local strings have an extra complication over global defects, which stems from 
the fact that we are unable to simulate the underlying field theory Instead, we 
approximate the true dynamics with line-like rclativistic strings using the Nambu 
equations of motion. We used a previously developed implementation ]|, |?J of this 
algorithm which simulates the network of strings neglecting damping effects due to 
the expansion of the Universe. The velocities of the strings segments in the network 
can be constrained to be integers due to the gauge conditions used to discretise the 
equations of motion and this enables the code to obtain very high accuracy with 
extremely fast computation. 

To simulate the extraction of energy from the system due to the decay of the 
string loops into gravitaional radiation, loops of a minimum size are excised from the 
simulation at each timestep. This also ensures that the network scales with respect 
to the conformal time r\ which enables us to extend the dynamical range of the 
resulting correlation functions beyond the limited range covered in the simulation. 

We performed simulations with box sizes ranging from 128 3 to 600 3 , with a cut- 
off on the loop size of two links. Realisation averages were carried out with 256 3 
boxes once it was determined that the general form of the correlators scaled very 
accurately with box size. To evaluate the UETCs from the simulations we selected 
times in the range 0.1 A~ < t < N/4, where N is the box size, when we were sure 
that the string network was scaling, and when boundary effects are still excluded by 
causality. From this we obtained the time evolution of the network's stress-energy 
tensor M „ at each point on the lattice using 

M „(x) =nj da{X^X» - A^JT)5 3 (x - X(<r, if)) (8) 

where X (17,77) is the spacetime coordinate of the string with a a parameter along 
the string. SVT decomposition was carried out on the 10 independent components 
of the Fourier transform of the stress-energy tensor at each timestep. This resulted 
in all the SVT components being computed directly from the simulation without 
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making assumptions on energy conservation and on the details of energy dissipation 
from the string network. The drawback of obtaining all the components directly 
in such a manner is that the process becomes computationally intensive for even 
modestly sized simulations e.g. 256 3 . By cross-correlating the decomposed stress- 
energy components from a central time with those from all the stored timesteps the 
14 independent UETC's f XY (k,rj,rj') were computed. 

Standard codes for structure formation [fto[ require the square root of coherent 
correlators for sourcing the inhomogcneous differential equations but this is not 
straightforward in general defect scenarios because the correlators are incoherent 
and therefore do not factorise. To overcome this, following [ffl], we expand the matrix 
of scaling functions onto a basis of coherent functions {krf} by diagonalising the 
matrix of correlators 

i 

where AW are eigenvalues. Each coherent, weighted eigenmode is then fed into 
the standard codes and the total Ci and -P(fc) are then the convergent sums of the 
resulting components. 



4. RESULTS 

The most striking difference between our work on local strings and previous 
work on global defects is that we find the Ooo term comes to dominate over all 
the other components and thus vector and tensor components contribute much less 
significantly in our results. The string anisotropic stresses are in the predicted |8| 
ratios |6 S | 2 : |6 y | 2 : |6 T | 2 of 3 : 2 : 4, as fcr -» 0, however |6 o| 2 > |6 S | 2 . Also we 
find that the peak in the energy density power spectra occurs at fcr ss 20 which is 
well inside the horizon. In general this leads to a higher Doppler peak than observed 
in SI- 

For a full treatment we must also include in the calculations the background fluid 
into which the network of strings dumps energy. One of the greatest unknowns in 
the dynamics of cosmic strings is the channel through which this energy is discarded 
and this is quantified by assuming an equation of state for the extra fluid, of the form 
p x = w x Q x Q , where w x is varied in the range < w x < 1/3 so as to account for 
the non-relativistic case through to the relativistic case. 

In Fig.l we plot y/[£(£ + l)Ce/2n], setting the Hubble constant to H a = 50 Km 
sec -1 Mpc -1 , the baryon fraction to fib = 0.05, and assuming a flat geometry, no 
cosmological constant, 3 massless neutrinos, standard recombination, and cold dark 
matter. The most interesting feature is the presence of a reasonably high Doppler 
peak at £ — 400 — 600, following a pronouncedly tilted large angle plateau. This 
feature sets local strings apart from global defects. It puts them in a better shape 
to face the current data. 

In Fig. 2 we plotted the CDM power spectrum P(fc) together with experimental 
points as in M. The normalization has been fixed by COBE data points. We 
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Fig.l The CMB power spectra predicted by cosmic strings decaying into loop and radiation 
fluids with W X = 1/3,0.1,0.01,0. We have plotted (1(1 + l)Ct/2n) 1/2 in flK, and 
superposed several experimental points, fluid. 
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Fig. 2 The power spectrum in CDM fluctuations for cosmic strings, with W J 
0.01,0.1, 1/3. We plotted also the standard CDM scenario prediction and points inferred 
by Peacock and Dodds from galaxy surveys. 
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see that the peak of the spectrum is always at smaller scales than standard CDM 
predictions, or observations. The CDM rms fluctuation in 8 /i _1 Mpc spheres is 
a$ = .42, .61, 1.8 for w x = 1/3,0.1,0.01. Hence relativistic decay products match 
well the observed as ~ 0.5. On the other hand in 100 /i _1 Mpc spheres one requires 
bias &ioo = ofQQ a /(7ioo = 4.9, 3.7, 1.6 to match observations. 

In this work we have calculated the CMB and LSS power spectra arising from 
local cosmic string scenarios using the largest string simulations to date. The spectra 
account for a variety of energy decay mechanisms and the result we want to highlight 
is the significant difference between structure formation processes arising from global 
and local defects which show that the 100 h~ x Mpc -1 bias problem and low Dopplcr 
peaks are not as generic effects as was previously thought. 
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